Mon. Not. R. Astron. Soc. OOO.rflfHlO Printed 29 March 2010 (MN KTbX style file v2.2) 



The transition from population III to population II-I star formation 



Umberto Maio 1 ' 2 * Benedetta Ciardi 1 , Klaus Dolag 1 , Luca Tornatore 3 , Sadegh Khochfar 2 

1 Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Strafie 1, D-85748 Garching b. MUnchen, Germany 

2 Max-Planck-Institut fur extraterrestrische Physik, Giessenbachstrafie 1, D-85748 Garching bei MUnchen, Germany 

3 Dipartimento di Astronomia dell'Universitd di Trieste, via Tiepolo 11, 1-34131 Trieste, Italy 



...(draft) 



ABSTRACT 

We present results from the first cosmological simulations which study the onset of pri- 
mordial, metal-free (population III), cosmic star formation and the transition to the present- 
day, metal-rich star formation (population II- 1), including molecular (H2, HD, etc.) evolution, 
tracing the injection of metals by supernova? into the surrounding intergalactic medium and 
following the change in the initial stellar mass function (IMF) according to the metallicity of 
the corresponding stellar population. Our investigation addresses the role of a wide variety of 
parameters (critical metallicity for the transition, IMF slope and range, SN/pair-instability SN 
metal yields, star formation threshold, resolution, etc.) on the metal-enrichment history and 
the associated transition in the star formation mode. All simulations present common trends. 
Metal enrichment is very patchy, with rare, unpolluted regions surviving at all redshifts, in- 
ducing the simultaneous presence of metal-free and metal-rich star formation regimes. As a 
result of the rapid pollution within high-density regions due to the first SN/pair-instability 
SN, local metallicity is quickly boosted above the critical metallicity for the transition. The 
population III regime lasts for a very short period during the first stages of star formation 
(~ 10 7 yr), and its average contribution to the total star formation rate density drops rapidly 
below - 1CT 3 - 1(T 2 . 
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1 INTRODUCTION 

The standard paradigm of cosmic struct ure formatio n relies on 
the classical approach of Jeans' theory jjeanslll902h applied to 
primordial matter fluctuations in the frame of an expanding Uni- 
verse. Cosmological models for struct ure formation have been de- 
veloped since several dec ades (e.g. iGunn & Got3ll972l ; |Peeblesl 
1 1974 IWhite & Reeil 19781) and the overall picture agrees with a 
"flat" Universe where "cold-dark matter" is the dominant frac- 
tion of matter and the "cosmological constant", A, is the domi- 
nant fraction of the cosmological energy density. Baryonic struc- 
tures arise from in-fall and condensation o f gas into dark-matter 
potential wells. In particular, it seems (e.g. iMaiolino et alj[2007h 
that molecular gas could account for a significant fraction of the 
dynamical mass of earl y objects. Recent det erminations of the cos- 
mological parameters (Ko matsu et al.l l2010) suggest a present-day 
expansion rate Ho — 70 km/s/Mpc (in units of 100 km/s/Mpc 
this parameter becomes h ~ 0.70), a total-matter density param- 
eter Ho™ ~ 0.272, with "baryonic" component Qot — 0.0456, 
and a cosmological-constant density parameter £Iqa — 0.728. 
The primordial power spectrum of perturbation is well fitted by 
a power law with index n ~ 0.96 and normalization via mass 
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variance within 8 Mpc//i radius as — 0.8. As a reference, it 
is common to define the standard ACDM model the one with 
the following parameters: Ho = 70 km/s/Mpc, fiom = 0.3, 
n b = 0.04, fi A = 0.7, n tot = 1.0, n = 1, erg = 0.9. 
In the frame of cosmic evolution, it is believed that structure for- 
mation takes place from the growth of primordial fluctuations in 
matter density. These would contract and collapse allowing gas 
cooling and the subsequent build-up of stars. The determination 
of the properties of early stars and their effects on the follow- 
ing baryonic-structure formation episode s is still a problem un - 
der debate (for a complete review see e.g. ICiardi & Ferrarall2005l) . 
Nevertheless, it is reasonably well established that the very first 
generation of stars should be characterized by massive objects 
with typical masses much larger than the presently observed ones 
dSchwarzschild & Spitzeill 19531: IWooslev & Weaver! 1 19951: ILarsoij 
ll998l:IChiosill2000l : lHeger & Wooslevll200ll2008h These primor- 
dial stars (population III stars) are formed out of a pristine envi- 
ronment, where the cooling agents are limited to primordial, H- 
based molecules only, i.e. H2 and HD, which are able to cool the 
gas down to temperatures of ~ 10 2 K. Therefore, the mass of 
primordial stars should be relatively large and their spectrum, bi- 
ased towards such large objects, is commonly ref erred to as "top- 
heavy" initial mass function (IMF) dLarson|[l998|) . These features 
imply very short lifetimes (up to ~ 10 6 years only) and final 
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death mostly into black holes dHeger & Wooslevll2002h . The only 
mass range where primordial stars can explode as pair-instability 
supernova s (PISN) and pollute th e surrounding medium is [140- 
260] M© dHeger & Wooslevi2002r) . Nevertheless, also at masses of 
~ 100 — 140 Mq there can be some mass loss before collapse 
(pulsational pair SN). The key uncertainty here is primary nitrogen 
production and the dredge up of carbon and oxygen. In particular, if 
the stellar atmosphere is highly CNO enhanced there may be sub- 
stantial mass loss, but zero metallicity should still be a good first 
approximation for such stellar flows (S. Woosley, private commu- 
nication). 

Despite the many uncertainties on their characteristics, population 
III (popIII) stars have an important impact on the evolution of the 
intergalactic medium (IGM), since they initiate the metal pollution 
of the IGM, with consequent change of its chemical composition 
and cooling properties (chemical feedback). Therefore, star forma- 
tion events in enriched regions will happen in completely different 
conditions, because metals allow further cooling and fragmenta- 
tion to smaller scales. This results in an initial stellar-mass function 
peaked at lower masses, similar to the nowadays observed Salpeter- 
like IMF i Salpet eill 19551) for population II-I (popII-I) stars. 
A very debated issue is the understanding of the transition from 
the primordial popIII star formation regime to the standard popll- 
I regime. In this respect, there are evidences for the existence of 
a critical metallicity, Zcrit, at which the modal i ties of star for- 
matio n allow such transition dBromm et al.ll20oH ; ISchneider et al.l 
2002): Z cr it is the metallicity at which the metal cooling function 
dominates over the molecular one. In this case, the fragmentation 
process becomes highly enhanced, but the exact value is not well- 
established, yet. Different studies sugg est discrepant values with 
Z cr it varying bet ween ^ 10 6 Zq (e.g . lSchneider et al .120061) and 
~ 10" 3 Z Q (e.g. lBromm & Loebll2003lf . 

The main uncertainty is the presence of dust at high redshift. To 
our knowledge, dust is produced and injected in the ISM mainly 
by the low-mass stars in the AGB phase. This would imply no 
dust production at very-high redshift, due to "long" stellar life- 
times: a 8 Mq star has a life of ~ 0.1 Gyr, comparable with the 
age of Universe at z ~ 15 — 20. None the less, the presence 
of large amounts of dust and heavy elements has been detected 
at moderately-high redshift, i .e. at z > 6, when the Universe is 
young er than ~ 1 Gyr (e.g. Bertoldi et al. 2003; Maiolin o et al.l 
2004). This suggests that dust production must have occurred pri- 
marily in the ejecta of supernova explosions, which are the final 
fate of massive, short-lived stars. If, indeed, supernovas could be 
able to induce dust production then star formation in the early 
Universe and the level of Z cr i t would be strongly influenced, as 
well. Enrichment by dust would not impact the thermal struc- 
ture of the IGM, but its presence would alter the whole star for- 
mation process. Indeed, when the metallicity of star forming re- 
gions is still below ~ 10 -6 Zq the only relevant coolants are 
molecules, mainly H2 and HD, while above ~ 10~ 4 Zq — 10 -3 Zq 
gas coolin g is fully domin ated by metal fine- structure transi- 
tions (e.g. [Maio et ai] |2007l) and cloud fragmentation can hap- 
pen down to sub-solar scales dBromm et al. 2001; Schneid er et al.l 
120021 12004) . In between, cooling capabilities depends mainly on 
the amount of metals depleted onto dust grains. The mass con- 
tent locked into dust grains, almost independently from the exact 



mass of SNII/PISN p rogenitor, reaches 2% — 5% of the parent 
mass for type-II SN dKozasa et all 1 199 it iTodini & Ferraral boOll : 
iNozawaet al.ll2003l;lBianchi & Schneidej|2007l). an d 15% - 30% 
for PISN dNozawa et alj2003l :ISchneider et al.l2004l) . in a few hun- 
dred days after the beginning of the explosion. Despite the huge 
difference between the dust mass fractions deriving from SNII and 
PISN, the dust-to-metal mass ratios (or "depletion factors") turn out 
to be ~ 0.3 — 0.7 in both cases and, for metal poor SNII progenitors 
with masses of 25 Mq — 30 Mq, it is up to ~ 1. The natural con- 
clusion is that first massive stars could spread in the early Universe 
a lot of dust during their final phases (as either SNII or PISN). 
In the present work, we aim at investigating the birth of the first 
stars, the following cosmic metal enrichment from their explosive 
death and the transition to the standard, presently observed, star 
formation regime. We focus on how such transition and the as- 
sociated features are affected by different choices of Z cr it, IMF, 
post-supernova metal yields, star formation density thresholds, box 
dimension and resolution. 

Throughout the paper we will refer to "population III regime" and 
"population II-I regime" when the metallicity is below or above the 
critical level, respectively. In addition, we will adopt the standard 
ACDM cosmological model. 

The paper is organized as follows: we describe the code and the 
simulations in Sect. [2] results are given is Sect. [3] and parameter 
dependence is addressed in Sect. l4l d4~Tl 14.21 and 14.31 ). We discuss 
and conclude in Sect. [5] 



2 SIMULATIONS 

The simulations were performed in the frame of the standard 
ACD M cosmological model, by using the Gadget-2 code tSpringel 
2005). The initial conditions for the reference simulations were 
generated sampling 0.7 Mpc//i, comoving (but see also discus- 
sion on larger boxes in Sect. 14.3b of the cosmic fluid at redshift 
2 = 100, with 320 3 particles per gas and dark matter species, hav- 
ing masses of 116 M0//1 and 755 Mq /h, respectively. We perform 
the numerical simulations including prim ordial molecular evolu- 
tion (Maio et al. 2006, 2007; Maio 2008), star formation prescrip- 
tions (M aio 2008; Maio et al. 2009), and metal pollution from stars 
( Torna tore et alj2007l) . We choose the stellar IMF according to the 
value of the gas metallicity Z: if Z < Z cr u a top-heavy IMF is 
assumed, otherwise a Salpeter one. 

More in detail, we include in the code the whole set of chemistry 
reactions leading to molecule creation or destruction and metal 
pollution, with the corresponding cooling functions. We follow 
the abundances of the H- and He-based chemistry, involving e~, 
H, H+, H~, He, H e+, He++; H 2 -chemistry following H 2 , H+ 
dYoshida et al]|2003l) and de uterium chemistry with D, D + , HD, 
and HeH + dMaio et all 20071) . While in principle an arbitrary num- 
ber of metal s can be accounted fo r by the scheme for metallic- 
ity tracking dTornatore et al.ll2007l) . in this work we focus on C, 
O, Mg, S, Si, and Fe, since they are the most abundant metals 
produced during stellar evolution. We assume, as initial chemical 
composition, a primordial neutra l gas with residual electron and 
H+ fractions dGalli & Pall ;]| 19981) x e - ~ z H+ ~ 4 • 10~ 4 , H 2 
fraction xu 2 = 10 -6 , H^ fraction x n + = 3 ■ 10 -21 , D frac- 
-5 n+ f™t; nn ™ — a . in— 9 



1 W e adopt Zq ~ 0.0201 dAnders & Grevess3ll989l ; lGrevesse & Saliva] 
1998). See, however, Asplund et al. 12009) for a recently updated value of 
Zq ~ 0.0134 



tion xt, = 3.5 • 10 , D + fraction x D + = 4 • 10 , HD fraction 
x'hd = 7- 10~ 10 , HeH+ fraction x HcH+ = 10~ 14 . Individual met- 
als are set to zero initial abundance. 

The natural gas evolution, with the initial in-fall into the dark- 
matter potential wells and the shock-heating phase, is fully self- 
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Table 1. Parameters adopted for the simulations. 



Model 


box side 


number of 


mean inter-particle 


M g as[M e /h] 


M dm [M Q /h] 


Z 


crit 


"H,th 


popIII IMF 


Massive 




[Mpc/A] 


particles 


separation [kpc/h] 










[cm" 3 ] 


range [Mq] 


yields 


Zcrit3 


0.7 


2 X 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10" 


3 z Q 


70 


[100, 500] 


HW02 1 


Zcrit4 


0.7 


2 X 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


w- 


4 z Q 


70 


[100, 500] 


HW02 1 


Zcrit5 


0.7 


2 x 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10- 


5 Z e 


70 


[100, 500] 


HW02 1 


Zcrit6 


0.7 


2 x 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10- 


6 z e 


70 


[100, 500] 


HW02 1 


Zcrit4-r8-100 


0.7 


2 x 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10" 


4 z Q 


70 


[8, 100] 


HW08 2 


Zcrit4-r0.1-100 


0.7 


2 x 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10- 


4 z & 


70 


[0.1, 100] 


HW08 2 


Zcrit4-HW8-40 


0.7 


2 x 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10" 


4 z Q 


70 


[0.1, 100] 


HW08 2 


Zcrit4-WW8-40 


0.7 


2 x 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10- 


4 z & 


70 


[0.1, 100] 


WW95 3 


Zcrit4-th7 


0.7 


2 x 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10" 


4 z Q 


7 


[100, 500] 


HW02 1 


Zcrit4-thl 


0.7 


2 x 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10- 


4 z Q 


1 


[100, 500] 


HW02 1 


Zcrit4-th01 


0.7 


2 x 320 3 


2.187 


1.16 x 10 2 


7.55 x 10 2 


10- 


4 z @ 


0.1 


[100, 500] 


HW02 1 


Zcrit4-5th7 


5.0 


2 x 320 3 


15.62 


4.23 x 10 4 


2.75 x 10 5 


io- 


4 z e 


7 


[100, 500] 


HW02 1 


Zcrit4-10th7 


10.0 


2 x 320 3 


31.25 


3.39 x 10 5 


2.20 x 10 6 


10- 


4 z Q 


7 


[100, 500] 


HW02 1 
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2 Data tables from 

3 Data tables from 
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Wooslev & Weaver Jl995h 



consistently followed through the attainment of the isothermal turn- 
over and the subsequent cooling regime. This stage is very im- 
portant at early times, as it represents the state when gas cool- 
ing balances and takes over heating. Resolving the isothermal 
phase in primordial-chemistry simulations is fundamental to ac- 
count for the early molecular cooling of the gas (as already dis- 
cussed and shown in iMaio etall 2009). For this reason, we assume 
that stars are formed roughly at the end of the gas cooling pro- 
cess, when a suitable physical density threshold for early st ar for- 
mation of ~ 10 2 cm" 3 (more exactly 70 cm" 3 ) is reached dMaiol 



120081; IMaio et alj|2009h. The often used low-density threshold of 
~ 0.1 cm -3 (e.g. iKatz et al. I ll996t) woul d induce spurio us time 
anticipation of the onset of star formation jMaio et all 2009). Here- 
after, stars evolve and die according to their characteristic lifetimes, 
possibly expelling newly created metals into the surro unding gas 
and m odifying its cooling properties e ither at low (e.g. Maio et all 
20071) and at high temperatures (e.g . ISutherland & Dopitdl 19931 : 



Wiersma et al 2009; Cantalu poll201 



For the very first bursts (and those happening in gas of metallic- 
ity smaller than Zcrit) we adopt a "Salpeter-like" IMF, in the mass 
range [100, 500] Mq, with a slope of —1.35. The stellar life-times 
lie between ~ 3 • 10 6 yr, for the 100 Mq stars, and ~ yr (instan- 
taneous death), for the hundreds-solar-masses stars. As mentioned 
in the Introduction, the relevant mass range for metal pollution is 
approximatively [140, 260] Mq, because of the formation of PISN. 
We point out that in our approach each star forming SPH particle 
can be assimilated to a simple stellar population, with coeval stars 
formed out from the same environment. Once the particle (i.e. stel- 
lar population) metallicity reaches the critical value, Zcrit, popII-I 
star formation is assumed to set in, in place of popIII, and subse- 
quent stars will have masses in the range [0.1, 100] Mq, distributed 
according to a Salpeter IMF. PopII-I stellar life-times are usually 
much longer than those of popIII, ranging between ~ 2 ■ 10 10 yr, 
for the 0.1 M stars, and ~ 3 ■ 10 6 yr, for the 100 M . 
Yields of popII I stars (see discussion in Sect. 14. 1 1 ) are assumed 
accordingly to iHeger & Wooslev! d2002l) . while for popII-I stars 



we follow Wooslev & Weaverl \ 19951) for massive stars (SNII), 
Ivan den Hoek & G roenewegen ( 1997) for ma s s loss of low- and 
intermediate-mass stars, and Thielema nn et alj J2003h for SNIa. 
In order to describe properly stellar evolution, it is important to take 
into account feedback effects which deal with those physical pro- 
cesses that are not resolved because of numerical limitations. In our 
simulations star forming particles are subject to wind kinetic feed- 
back with velocities of 500 km/s. These particles produce metals 
which are smoothed over the neighbours according to the SPH ker- 
nel. 

We perform four numerical simulations considering critical metal- 
licities of Z crit = 1O" 3 Z , W~ 4 Z Q , W~ 5 Z Q , W~ 6 Z Q , which 
span the full range found in the literature (see Introduction). In 
Sect. 14. 1114.21 and !4.3l we will also investigate how the results are af- 
fected by a different choice for the popIII IMF, the post-supernova 
yields, the critical star formation threshold and the box dimen- 
sion/resolution. A schematic description of the simulations is given 
in Table[0 



3 RESULTS 

In this section, we present results related to the early phases of cos- 
mological metal enrichment and the transition from popIII to popll- 
I regime. The results refer to the 0.7Mpc//i side box described in 
the previous section, with varying Z cr it, the most significant pa- 
rameter in the calculation. 



3.1 Metallicity distribution 

In order to have a visual representation of the simulations we show 
maps at redshift z = 11 for the projected metallicity distribution 



2 For a different metal/feedb ack approach, se e also treatments based on 
thermal feedback (e.g. Scanna pieco et alj2005t) . 



4 



Umberto Maio et al. 




-3.0 



■ -4.0 , 



-5.0 



-3.0 



-4.0 



Figure 1. Metallicity maps at redshift z = 11 for the 0.7Mpc//i side box, with Z crit = 10 z Zq (top-left panel), 10 
(bottom-left panel), and 10 — 6 Zq (bottom-right panel). The maps are obtained by projecting the whole simulation box. 



1 Zq (top-right panel), lO -5 ^ 



(Fig-ID f° r al l trle cases considered: Z cr i t = 10~ 3 Zq (upper-left 
panel), lO _4 Z0(upper-right panel), lO _5 i?0(bottom-left panel) 
and 10~ 6 Zq (bottom-right panel). The metallicity is defined as the 
metal mass fraction of each SPH particle: Z = mz/nitot, where 
mz and mtot are respectively the mass of metals and the total mass 
of each SPH particle. 

Metal enrichment is very patchy and there are regions which are 
very strongly polluted and others which preserve their pristine, 
metal-free composition. The main pollution events happen around 
the densest regions, where star formation takes place. We stress that 
the differences in the distribution and level of metallicities are only 
due to the value of Z cr it adopted. In particular, the metallicity of 
the gas reaches a slightly higher value and wider distribution for 
the Z cr it = lO~ 3 /70 case than for the other cases. This is due to 
the longer time needed to reach the critical level, Z CT it- Once this 
happens, popII-I star formation sets in and further delays metal pol- 
lution, because the SNII from popII-I stars release about one fourth 
of metals per unit mass than the popIII ones, on a timescale at least 
~ 10 times longer. When Z cr i t = 1O _6 Z0, pollution above the 
critical metallicity threshold is immediate and most of the regions 
rapidly move into the popII-I star formation regime. This transition 
gets delayed as Z cr it increases. 

In Fig. [2] we compare metal enrichment at different redshifts. The 
maps refer to z — 11 (upper row), 12 (central row) and 13 (bottom 
row), for Z cr it = W~ [i ZQ (left column) and W~ 5 Zq (right col- 
umn). The process is quite fast as inferred from the rapid growth of 
metal-enriched bubbles, and allows for high metal pollution within 



very short time (see also Sect. |3.2| l. 

To highlight the differences between popIII and popII-I star form- 
ing regions, we plot in Fig. [5] the corresponding two-population 
map!|j projecting the whole box. They show very well how the star 
formation regime evolves according to metallicity (see further dis- 
cussion in Sect. l3.3t . 

From a more quantitative analysis of the Figure, some distinctive 
behaviour emerges. For Z CT it = 10~ 6 Zq, almost all the polluted 
areas reach the critical metallicity already in ~ 0.1 Gyr (~ 98% 
at z = 13, ~ 99% at z = 12, and ~ 100% at z = 11). 
In the other cases we observe a similar quick enrichment: for 
Zcrit = W~ 5 Zq, 10~ 4 Zq and 10~ 3 Zq the percentage of areas 
with Z > Zcrit is respectively ~ 72%, 49% and 32% at z = 13; 
~ 86%, 70% and 39% at z = 12; and ~ 99%, 97% and 97% 
at z = 11. We note as well that, at fixed redshift, the fraction of 
popII-I star formation regions decreases with increasing Z cr it, be- 
cause of the higher enrichment level needed to allow the transition 
from Z < Zcrit to Z > Z crit . 

As the enrichment process is very localized, also the transition to a 
different star formation regime strongly depends on location. This 
means that it is possible to find coexistence of population III and 
population H-I, with stronger pollution in the central regions of star 
formation sites and weaker pollution in the outermost ones (see 
maps in Fig. Q]and in Fig. [2j, It is in fact evident that the pref- 



In this case, data are not smoothed on the edges of the pixels. 
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Figure 2. Comparison of the metallicity maps at redshift z = 11 (upper row), 12 (central row), and 13 (bottom row), for the 0.7Mpc/Ji side box with 
Zcrit = W~ 3 Zq (left column) and 10 -5 Zq (right column). The maps are obtained by projecting the whole simulation box. 



erential sites of popIII star formation are isolated areas and often 
peripheric regions not strongly contaminated by nearby bursts. Our 
findin gs are consis t ent wi th the study of the Z CT it = 10~ 4 Z@ case 
bv lTornatore et al.l d2007t) . We note that changes in Z cr it can only 
slightly alter the whole enrichment scenario, as just few regions 
which undergo popII-I star formation for Z aT n — W Z@ are un- 
dergoing popIII star formation for Z cr u = lO -3 /?© (e.g. Fig. [3j(. 
In fact, accordingly with the previous discussion, the star forming 
environment is chemically "sterilized" by one or a few subsequent 
popIII generations. 



3.2 Metallicity evolution 

From the previous results, it emerges that metal enrichment is a 
very patchy and inhomogeneous phenomenon. Never the less, as 
an indicative estimate of the global enrichment history of the Uni- 
verse, we present the evolution of the average metallicities found 
in the simulations (Fig.[4j. Following previous discussions though, 
we expect large fluctuation jf] around the mean values reported here. 
In particular, for each Z cr i t (labeled on the upper-right corner and 



4 Variances are comparable to the averages. 



6 



Umberto Maio et al. 



2 cn - t =10-% un 



z=1 1 .00 



1 1 .00 





Figure 3. Compai'ison of the two-population metallicity maps at redshift z = 11 (upper row), 12 (central row), and 13 (bottom row), for the 0.7 Mpc/fe side 
box with Z cr it = 10~ 3 Zq (left column) and 10 — 5 Zq (right column). The yellow areas refer to polluted regions with metallicity higher than Z cr n, while 
the red ones to those with metallicity lower than Z ar u. The maps are obtained projecting the whole box. 



shown by the horizontal dot-dot-dot-dashed line in each panel), 
we show the maximum metallicity reached at any redshift (dotted 
lines), the average metallicity of the enriched regions (dot-dashed 
lines), the average metallicity in the whole simulation box includ- 
ing primordial unpolluted zones (solid line) and the corresponding 



individual metallicities for the different metals tracked by the code 
(dashed lines). 

The first metals spread during the final stages of stellar evolu- 
tion leads to typical metallicities in the surrounding ISM of Z ~ 
(10 -5 — 10~ 4 ) Zq, which, shortly after, rise to values of Z ~ 
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Figure 4. Metal evolution as function of redshift for the 0.7Mpc/h side box simulations, with Z cr it = 10 _3 Zq (upper-left panel), 10 — 4 Zq (upper-right 
panel), 10 — 5 Zq (bottom-left panel), and 10 _6 Zq (bottom-right panel) case. The magenta horizontal dot-dot-dot-dashed line indicates, in each panel, the 
critical metallicity; the dotted line is the maximum metallicity; the dot-dashed line is the average metallicity of the spread metals; the solid line is the total 
metallicity averaged over the whole simulation box and the dashed lines the corresponding individual metallicities: oxygen (green), carbon (blue), magnesium 
(red), sulphur (purple), silicon (pink), iron (yellow), other metals (cyan) - see also text in Sect. [5] 



(1CT 4 - 10~ 3 ) Z Q , between redshift z ~ 16 and z ~ 15. There- 
fore, the critical metallicity Z cr i t is easily overtaken, despite its 
precise, actual determination. The leading element is always oxy- 
gen, as it is the most abundant one produced by SNII and PISN 
explosions. 

For higher Z cr u the popIII regime lasts longer, so that more mas- 
sive star explosions take place which can pollute the medium up 
to higher metallicities, before popII-I star formation regime sets in. 
This results in a quicker Z increase in the very early stages (the dif- 
ference between the Z CT it = Zq and Z cr i t = W~ 6 Zq cases 
is almost two orders of magnitude at z ~ 15). 
Broadly speaking, a metallicity of 10~ 6 — 10~ 5 Zq is locally 
reached after only ~ 10 5 yr. In ~ 5 • 10 6 yr, it is possible to get 
Z ~ 10 -4 Zq or even Z ~ 10 -3 Zq. An average pollution of 
Z ~ Z C rit is always obtained in a few 10 7 yr (i.e. Az ~ 1, at 
z ~ 16) after the onset of star formation (Fig. [4}- Once the critical 
metallicity is reached, popII-I star formation sets in and contributes 
(on longer timescales) to metal enrichment, as well. 



3.3 Star formation from population III and population II-I 
stellar generations 

In Fig. [5] we plot the averag^f] popIII star formation rate, SFRm, 
and the average total star formation rate, SFR to t- For all the simu- 
lations the onset of star formation is at z ~ 16.3, when the Universe 
is about 2.3 • 10 8 yr old (independently from Z cr i t , which plays a 
role only in the following epochs). As shown in Fig. [5] at the very 
beginning popIII contribution is dominant, but rapidly fades away 
and the popII-I regime is established. 

For the Z cr it = 10 _3 Zq case, the contribution is higher and it de- 
creases with Z cr it. In fact, the higher Zcrit, the higher the number 
of star forming regions with Z < Zcrit ■ The total star formation 
rate density is only mildly influenced by the exact value adopted for 
Zcrit, meaning that popIII star formation does not affect the global 
behaviour significantly and the bulk of star formation is mainly led 
by popII-I stars. The reason for this is simply understood in terms of 
timescales, as popIII stars have shorter lifetimes (at most ~ 10 6 yr) 

5 The average is the arithmetic average of the star formation rates taken at 
each time-step and binned over 50 bins in redshift. 




Figure 5. Upper panel: redshift evolution of the average star formation 
rate densities (in Mqvt - 1 Mpc — 3 ) for the 0.7Mpc/ft side box simula- 
tions. The four upper lines are the total star formation rate densities for 
Z crit = 1CT 3 Zq (solid line), 1O" 4 Z (dotted line), lCr 5 Z Q (dashed 
line), and 10 — 6 Zq (dotted-dashed line). The four lower lines are the corre- 
sponding popIII contributions. Bottom panel: evolution of the ratio between 
the average popIII star formation rate density (SFRni) and the total star 
formation rate density (SFR tot ), SFRm/SFRtot- The bullets correspond 
to the redshift at which the average metallicity of the enriched particles 
equals the value of Z ar it (as in Fig. [4). 



Figure 6. Upper panel: redshift evolution of the average star formation rate 
densities (in M yr — 1 Mpc — 3 ) for the 0.7Mpc/ft side box simulations 
with Z cr a = 1O~ 4 Z . The upper lines are the total star formation rate 
densities for different popIII IMF mass range: [0.1 Mq, 100 Mg] (solid 
lines), [8 M ,1OO Mg] (dashed lines), and [100 M ,5OO M Q ] (dot- 
ted lines). The lower lines are the corresponding popIII contributions. Bot- 
tom panel: evolution of the ratio between the average popIII star formation 
rate density (SFRni) and the total star formation rate density (SFR to t), 
SFRm/SFRtot- The bullets correspond to the redshift at which the aver- 
age metallicity of the enriched particles equals the value of Z cr n. 



and heavily pollute the medium up to Z crit via PISN explosions. 
Therefore, after the first bursts, it is much easier to match the con- 
ditions for standard popII-I star formation, rather than metal free 
star formation. The period over which star formation is dominated 
by popIII is slightly longer in the Z C rtt = 1O _3 Z case and de- 
creases gradually with Z crit . For Z crit = 1O" 3 Z , ~ 100% of the 
total SFR is due to popIII stars in the redshift interval z ~ 15 — 16 
(equivalent to a time interval of ~ 2 ■ 10 7 yr), while the percent- 
age decreases to an average value of roughly 0.1% at z ~ 12. In 
the Z cr u = 1O _6 Z case, there is a sudden drop below ~ 10% 
already at z ~ 15 — 16 and down to ~ 0.01% at z ~ 12. The 
intermediate cases are bracketed by the former two regimes. These 
trends are consistent, because the time to pollute the IGM to lower 
Z cr it is shorter. The resulting behaviour - the smaller Zcrit, the 
earlier the transition from popIII to popII-I dominated star forma- 
tion - is expected, because, as mentioned, the time needed to pol- 
lute the IGM up to lower Z CT it is shorter. 

The dots on the lines of the average fractions of popIII star forma- 
tion rate densities, point to the redshifts when the average enrich- 
ment reaches the level of Z cr it (i.e. the abscissa of the intersecting 
point between the dotted line and the horizontal line in each panel 
of Fig.[4]l for the different cases. At those redshifts, the average con- 
tribution of popIII star formation has already dropped of more than 
one order of magnitude. For the Z cr it = 10~ 6 Zq case, Z cr it is 
reached in the polluted regions when SFRm/SFRtot ~ 1 and this 
ratio steeply decreases afterwards. In the cases with higher Z cr it, 
the popII-I regime develops simultaneously with the popIII regime, 
and the average popIII contribution to the total star formation rate 
when Z ~ Z cr it is below ~ 10%. 



4 PARAMETER DEPENDENCE 

In the following, we will investigate the effects on our results of 
the main parameters (other than Zcrit) adopted in the simulations 
and choose as reference the one with box dimension 0.7Mpc/ft, 
Zcrit = 10~ 4 Zq, nn,th = 70 cm -3 and Salpeter slope for popIII 
stars in the range [1OOM , 5OOM ]. We have also checked the ef- 
fects of adopting a different imple mentation of SPH (the relative 
pressure SPH, rpSPH, according to lAbej|201cb . which only con- 
siders pressure gradients in force calculations and performes closer 
to grid-based codes. We found no strong deviations in the g eneral 
gas st ate (as also expected from the discussion in Sect. 3.8 of Abel 
2010) and even in the star formation rates. Metal pollution at early 
times is slightly more scattered, due to the capability of the rpSPH 
implementation of better resolving instabilities. Overall, there are 
no substantial changes in our results. 



4.1 Different population III IMF's and yields 

In order to account for the uncertainties on the popIII IMF, we ex- 
plore first the role of the slope of the popIII IMF, which so far has 
been fixed to the value s = —1.35. 

There are no robust arguments (neither theoretical nor observa- 
tional) supporting the above assumption and different slopes or 
shapes are not excluded. None the less, by changing slope, the ex- 
pected fraction of PISN is not dramatically modified, as its value is 
always around ~ 0.4 for any s between roughly —1 and —3. Thus, 
the results on metallicity evolution and star formation rates are not 
altered at all. 

Differently, the chosen mass range and yields of popIII stars could 
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Figure 7. Upper panel: redshift evolution of the average star formation 
rate densities (in M yr — 1 Mpc — 3 ) for the 0.7Mpc/ft side box simula- 
tions with Z cr n = 10 — 4 Zq. The upper lines are the total star forma- 
tion rate densities for different popIII metal yields, and fixed the popIII 
IMF mass range at [0.1 Mq . 100 M ]: primordial SN yields for stars of 
[8 Mq , 40 M ] (case Z=0 in lWooslev & Weavejl99^) (dashed lines! pri- 
mordial SN yields for stars of [8 Mm .40 Mml faeger & Wooslevll2008t) 
(dot-dashed lines), and p rimordial SN yields for stars of [8 Mq , 100 Mq] 
I Heger & Woosley 2008) (solid lines). The lower lines are the correspond- 
ing popIII contributions. Bottom panel: evolution of the ratio between the 
average popIII star formation rate density (SFRm) and the total star for- 
mation rate density (SFRtot), SFRni/SFRtot- The bullets correspond to 
the redshift at which the average metallicity of the enriched particles equals 
the value of Z cr a . 



short lifetime of stars with masses in the range [100 Mq, 500 Mq] 
(and in particular in the PISN range), which pollute the surrounding 
medium in a time-lag much shorter than the one needed to longer- 
living ~ 10 Mq stars. Therefore, metal pollution and the transition 
to popII-I regime happen earlier, with a residual popIII contribution 
of a factor of 3 or 4 smaller. The difference between the extreme 
case of stellar masses in the range [0.1 Mq, 100 Mq] and the case 
of stellar masses in the range [8 Mq, 100 Mq] is simply due to the 
different normalization of the IMF, since yields and massive-SN 
lifetimes are the same. Indeed, the SN fractions are ~ 0.14 and 1, 
respectively. Thus, as in the latter case 100% of the stars dies as 
SN, their enrichment results slightly higher, and the popIII contri- 
bution drops faster than in the former one (see also the bullets in 
the lower panel). 

Finally, we investigate the effe c ts of differ ent metal yields, 
accord ing to IWooslev & Weaveil d 19951) and [Heger & Wooslevl 
(2008). As the massive-SN mass range differs in the two works, 
in order to have a fair comparison between the yields only 
and to exclude normalization effects, we fix the popIII IMF 
range as [0.1 Mq,100 Mq] and the massive-SN mass range 
as [8 Mq,40 Mq ], In F ig. [7] we compare the results for the 



Woosley & Weaveil Jl995l) 's yields (dashed lines) and for the 
Heger & Woosley ( I2008l) 's ones (dot-dashed lines). For sake of 



completeness, we re-plot also the results for the IMF mass range of 
[0.1 Mq, 100 Mq] and t he massive-SN range of [8 Mq, 100 M ] 
(IHeger & Wooslevil2008h . Due to the longer lifetimes of 10 M - 
40 Mq stars (dashed lines and dot-dashed lines), with respect to 
PISN and 100 Mq star lifetimes (solid lines), the popIII regime 
can now last longer (compare bullets in Fig.|7]with the ones in Fig. 
[5]or in Fig.[6]l, and contribute up to 10% the total SFR at z ~ 11: 
some order of magnitude more than all the other cases. By directly 
comparing the effects of adopting different yields, in the same mass 
range, it comes out that their impact is negligible and there are no 
significant changes in the star formation history. 



play a more important role, since they directly affect metal pollu- 
tion and star formation history. The existence of popIII stars with 

mass es below ~ 10 2 Mq is supported by numerical simulations 

(e.g lYoshid3 120061: lYoshida et ail 120071 ; ICampbell & Lattanziol 



120081 ; ISuda & Fuiimotol I2010I) . which find that it is possible to 
fragment primordia l clou ds down to ~ 10 Mq. According to 
IWooslev & Weaveil j 1995b . Z = massive stars can die as SN 
if their mass is smaller than 40 Mq; above this limit, they col- 
lapse into black holes. Recent updates, though, suggest that such 
stars can die as SN even if their masses are larger, up to 100 Mq 
(Heger & Woosley 2008). We have investigated the implications of 
the uncertainties on our results. 

In Fig. [6] together with the results for our reference run (dot- 
ted lines), we show the trends for different popIII IMF mass 
ranges (while we keep the same slope), assuming metal yields 
for primordial, massive SN from IHeger & Wooslevl d2008ln . We 
plot star formation rates and the corresponding popIII contribution 
for popIII IMF mass ranges of: [0.1 M Q , 100 M ] (solid lines), 
[8 M Q , 100 M ] (dashed lines), and [100 M©,500 M ] (dotted 
lines). The latter case is our reference run and its popIII contribu- 
tion is always lower than the other ones. The reason is the extremely 



6 We use the tabulated post-supernova yields (Table 8 in Heg er & Wooslevl 
120081) for the 10 M - 100 Mq stars, explosion energy of 1.2 Bethe, 
standard mixing prescriptions, S=4 piston parameter (S. Woosley, private 
communication). 



4.2 Density threshold for star formation 

We now address how the transition from popIII to popII-I depends 
on the choice of the density threshold for star formation. Indeed, it 
is usual to assume that star formation takes place only in converg- 
ing flows where a given number density is reached. Because dif- 
ferent thresholds lead to different onsets of star formation (see also 
iMaio et alj2009l) . this could also affect the properties of the transi- 
tion. For this reason, we re-run the reference simulation (nn.th = 
70 cm -3 ) assuming density thresholds of nn,th = 7, 1, 0.1 cm -3 . 
Also in these cases the enrichment is very fast and the critical level 
in polluted environments is always reached between redshifts ~ 15 
and 18, very close to the value of our fiducial case, ~ 15 (see Fig. 
|4). Obviously, the higher the density threshold the later the onset of 
star formation, as the gas needs more time to condense. 
In Fig.[8]we show the star formation rates and the contribution from 
popIII stars for different density thresholds. The asymptotic am- 
plitude of the total SFR is not affected, but the contribution from 
popIII stars has different slopes for different nH.th: the higher the 
threshold, the steeper the evolution of SFRni/SFR to t, with this 
ratio decreasing of ~ 3 orders of magnitude in redshift intervals of 
Az ~ 15, 8, 6, 4, for n H ,tr = 0.1, 1, 7, 70, respectively. 
It is evident that the difference between the 70 cm" 3 , 7 cm -3 , and 
1 cm -3 cases is quite small, with star formation rates almost over- 
lapping since the beginning of the process. The reason is that such 
values fall in correspondence of the gas cooling branch, which 
is extremely fast and therefore the effects on the final results are 
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Figure 8. Upper panel: redshift evolution of the average star formation 
rate densities (in M0yr — 1 Mpc — 3 ) for the 0.7Mpc/h side box simu- 
lations with Z cr a = 10 — 4 Zq. The upper lines are the total star for- 
mation rate densities for different H-density thresholds for star formation, 
n n,th = 0.1cm -3 (dotted lines), lcm" 3 (short-dashed lines), 7cm -3 
(long-dashed lines), and 70 cm -3 (solid lines). The four lower lines are 
the corresponding popIII contributions. Bottom panel: evolution of the ratio 
between the average popIII star formation rate density (SFRni) and the 
total star formation rate density (SFR to t), SFRm/SFRtot- The bullets 
correspond to the redshift at which the average metallicity of the enriched 
particles equals the value of Z cr j t . 

very little. For the 0.1 cm -3 case, instead, we find discrepancies at 
higher redshift due to the unresolved isothermal peak. This is due 
to the fact that (unphysical) low-density threshold models convert 
gas into stars quite regularly and smoothly, when it is still shock- 
heating in the hosting dark-matter potential well. But, if high- 
density thresholds are used, the dense gas is converted into stars 
later on, only after the loitering, isothermal phase (during which 
it just gets accumulated at temperatures around ~ 10 4 K) and the 
catastrophic run-away cooling. So, large amounts of particles un- 
dergo star formation in a shorter time, after which the SFR catches 
up and converges to the low-density-threshold cases. In the mean- 
time, the critical metallicity is more quickly reached, as well. This 
also explains the differences among the slopes of the decrement. 
The higher the density thresholds, the higher the gas accumulated 
and catastrophically converted into stars, the higher the enrichment, 
and the steeper the ratio SFRm/SFRtot- 

The popIII/popII-I transition (see bullets) happens in all the cases 
when the contribution of SFRm has dropped roughly between 
10" 2 - 10~\ and at redshift z ~ 10 it is ~ 10" 3 . 
These conclusions support the fact that the transition from popIII to 
popII-I star formation regime is very fast, independently from the 
detailed numerical parameters, as long as the gas isothermal phase 
is resolved. 

4.3 Box size 

In this section, we discuss the effect on our results of different 
choices for the box size. In the next, we will show metallicity prop- 
erties (Sect. 14.3. It and star formation rate densities (Sect. I4.3.2t 



Figure 9. Average metallicities in the simulation of the box with a comov- 
ing side length of 5 Mpc//i and Z cr n = 10~ 4 Zq. The density threshold is 
7 cm -3 . The magenta horizontal dot-dot-dot-dashed line indicates the crit- 
ical metallicity; the dotted line is the maximum metallicity; the dot-dashed 
line is the average metallicity of the spread metals; the solid line is the total 
metallicity averaged over the whole simulation box; and the dashed lines 
the corresponding individual metallicities (see also Fig. [4). 



for simulations with the same parameters used for our reference 
simulation but with a box size in comoving units of 5 Mpc/h and 
10Mpc//i. Properties and resolution are also summarized in Ta- 
ble Q] Because of the lower resolution of such boxes (gas par- 
ticle masses are of the order ~ 10 4 — 10 5 M Q //i), we cannot 
adopt a very-high-density threshold (of ~ 10 2 cm -3 or more), 
therefore, to make a meaningful comparison we show results for 
n-H,th = 7 cm -3 . 



4.3.1 Metallicity 

The global behaviour is patchy and similar to the maps already 
shown for the 0.7 Mpc/h cases (Figs. Q] [2] [3). The only differ- 
ence is the larger sampling of star forming regions with a wider 
statistical significance. 

In Figs.[9]and[T0] we plot metallicity evolution for the simulations 
with 5 Mpc//i and 10 Mpc//i box side, respectively. In both cases, 
the critical metallicity is reached quite quickly and metal-rich re- 
gions have Z > Z cr u already by redshift ~ 10. The rapid at- 
tainment of Z cr it is independent of the resolution, which instead 
influences the onset of star formation. That is the reason why, in 
the 10 Mpc/h box side, star formation is slightly delayed. As in 
Fig-El we a ' so plot the maximum (dotted lines) and average metal- 
licity in the whole simulation box (solid lines) with the different 
metal contributions (dashed lines). We find that the first can eas- 
ily reach solar values, because of strong pollution due to the high 
PISN metal yields; the second has Z < Z cr i t at any redshift above 
~ 7. These fluctuations are reflected in the high patchiness of early 
metal spreading. The metal pollution histories do not show large 
discrepancies and, in spite of differences in resolution of more than 
one order of magnitude for the SPH masses, an overall agreement 
between larger (5 or 10 Mpc/h) and smaller (0.7 Mpc/h) boxes is 
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Figure 10. Average metallicities in the simulation of the box with a comov- 
ing side length of 10 Mpc/// and Z cr n = 1Q- 4 Zq. The density thresh- 
old is 7cm -3 . The magenta horizontal dot-dot-dot-dashed line indicates 
the critical metallicity; the dotted line is the maximum metallicity; the dot- 
dashed line is the average metallicity of the spread metals; the solid line is 
the total metallicity averaged over the whole simulation box and the dashed 
lines the corresponding individual metallicities (see also Fig. [4). 



reached. Thus, the general results are independent from the partic- 
ular cosmological sampling. 



4.3.2 Star formation from population III and population II-I 
stellar generations 

We will now discuss the behaviour of the star formation rates for 
the simulations previously introduced. 

The average star formation rate densities are shown in the up- 
per panel of Fig.[TT] for the boxes of 0.7Mpc//i (dashed lines), 
5Mpc//i (dotted lines), and 10Mpc//i (solid lines). The upper 
lines refer to the total star formation rate densities and the ones be- 
low to the corr esponding popIII s tar formation rate densities (data 
points are by iReddv e t al. 2008). Overall, they seem to connect 
fairly well with the observed trend at lower redshift. 
In the bottom panel, we present the ratios of popIII versus total star 
formation rate. Independently from the box size, the popIII contri- 
bution is unity at the very beginning (i.e. ~ 100% of the SFR is 
due to popIII) and, for all the cases, rapidly decreases below ~ 1% 
at redshift z ~ 11, and below ~ 0.1% at redshift z < 11: a time- 
lag much shorter than a Gyr. This means, consistently with what 
said in Sect. [3] that the bulk of popIII episodes happpens in a short 
time after the onset of star formation and it dominates the cosmo- 
logical SFR for a limited period at high redshift. At lower redshift, 
only isolated regions could still have pristine environment and host 
popIII star formation (which is none the less several order of mag- 
nitude smaller that the total SFR). 

A general feature of the large box simulations is the slight delay of 
the onset of star formation: this is due to the limited capabilities of 
the bigger SPH particles to condense, cool, and collapse. Indeed, 
molecular cooling is not efficient when the SPH masses are so high 
(~ 10 4 - 10 5 M o //i) and the primordial gas in-fall is not well re- 



Figure 11. Upper panel: star formation rate densities for the simulations 
with a comoving box side of 0.7 Mpc/ h (dashed lines), 5 Mpc//i (dotted 
lines) and 10Mpc//i (solid lines). The density threshold is 7cm -3 . For 
each of them, we show the total star formation rate density (upp er lines) 
and the popIII one (lower lines). Data points jReddv et alj |2008) are dif- 
ferent determinations of the star formation rate density: square symbols are 
the UV determinations, triangular symbols are the IR determinations. Bot- 
tom panel: evolution of the ratio between the average popIII star formation 
rate density (SFRni) and the total star formation rate density (SFRtot), 
SFRm/SFRtot. The bullets correspond to the redshift at which the aver- 
age metallicity of the enriched particles equals the value of Z cr n ■ 

solved. The qualitative agreement of these results is encouraging, 
though. 



5 DISCUSSION AND CONCLUSIONS 

In this work, we perform numerical simulations of early struc- 
ture formation, including both primordial, molecular evolution and 
metal enrichment from stellar death, to study the transition from 
an early, massive star formation mode to a more standa rd one, reg- 
ulated by the gas metallicit y. We follow (see details in lMaio et al.l 
l2007l ; lTornatore et al.l2007h the abundances of e~ , H, H+ , H" , He, 
He+, He++, H 2 , H+, D, D+, HD, HeH+, C, O, Mg, S, Si, Fe, and 
use different initial stellar-mass functions, metal yields, and critical 
metallicities for the transition from a popIII to a popll star forma- 
tion mode, Z cr it. At the present, our work is the only one dealing 
with detailed chemical evolution, from early molecule creation to 
the later stages of star formation, and, at the same time, allowing 
to trace simultaneous, different stellar populations, according to the 
underlying metallicity. 

Indeed, early structure formation can be accurately modeled only 
with a proper treatment of both the chemistry of primordial 
molecules and the metal enrichment. In fact, the main influence 
of chemical evolution on the following generations of structures 
is via metal pollution (chemical feedback). This event can com- 
pletely alter the cooling properties of the gas and thus the modal- 
ities of star formation, inducing a transition from a top-heavy to a 
standard IMF. The transition is believed to happen when the gas is 
enriched above Z cr it, which allows fragmentation below the typi- 
cal scales determined by primordial molecular cooling. Because of 



12 



Umberto Maio et al. 



our ignorance of the features of early dust formation and the lack 
of precision of many atomic and molecular data, the exact determi- 
nation of Z C rit is still elusive; reasonable values should range be- 
tween 10" 6 Z Q and 10~ 3 Z e dBromm et alj200ll ; lschneider et all 
120021) . 

In our simulations we assume four different values of Z cr i t in 
the above range: 10~ 6 , 10" 5 , 10" 4 , and 10~ 3 Z s . From our in- 
vestigation, some common features emerge. In general, metal pol- 
lution, independently from the parameters adopted in the simu- 
lations, is very patchy, with excursions of orders of magnitudes 
at all redshifts. T his is consistent with any simulation including 
metal enrichment |Raiteri et all 19 96; Gnedin 1998; M osconi et all 
200ll;lLia et al.l2002allbl:lKawata & Gibsoij2003l;lKobavashill2004l; 
Ricotti & Ostriker 2004; Scanna pieco et al.l 120051 ; iTornatore et al.l 
20071 : lOppenheimer et all l2009 ) . In addition, the gas is easily en- 



riched above Z cr u . For this reason, the average contribution from 
pristine, metal-free (or Z < Z cr it) stars to the total cosmic star for- 
mation density is dominant only in the very early phases of struc- 
ture formation, while it drops below ~ 10 -3 quite rapidly, after 
the explosion of the first pair-instability supernova? and their metal 
ejection. In fact, PISN explosions which follow the death of the 
first, metal-free or very-metal-poor stars, are the main responsibles 
for enriching the surrounding medium up to a minimum level of 
- 10" 4 Z Q . This means that nearby star forming regions have a 
very high probability of being polluted above Z cr it (as seen also in 
Fig-EJl, while popIII star formation can still occur farther away, in 
rare, isolated regions with pristine or low-metallicity gas (e.g. Fig. 
0. None the less, this would not be the dominant star formation 
regime. 

Our findings hold regardless of the value of Z cr u, the popIII IMF 
adopted, and the numerical parameters involved in modeling star 
formation (as shown in Sect. 14.11 14.21 and 14.3b . Differences are 
found if different IMF mass ranges for primordial stars are used. 
Because of the short life and the high metal yields of early, massive 
SN, the popIII regime contribute s, in any case, on ly slightly to the 
global SFR (as expected by e.g. iRicotti & Ostriker 2004), since the 
early pollution events quickly raise Z above Z cr it, independently 
from the detailed prescriptions. 

The simulations were performed using a standard ACDM cosmol- 
ogy, but slightly different parameters would not change the gen- 
eral picture. We have checked this, by running simulations with 
Clom = 0.26, n A = 0.74, Q pt = 0.0441, h = 0.7 2, ct 8 = 0.796, 
and n = 0.96 (WMAP5 data, Iffinshaw et al]|2008l) and found that 
the same results hold, albeit shifted by a Az ~ of a few to lower 
redshift (because of the smaller fiom). 

One last comment about the role of dust production from SNII 
and/or PISN: one of the main uncertainties in determining Z cr u. 
According to our findings, given the strength and rapidity of metal 
enrichment (see e.g. Figs. l4l [91 andl lOll. details about dust and its 
impact on Z cr u are not relevant for the general process and for our 
understanding of the transition from primordial, popIII regime to 
present-day-like, popII-I regime (even if they could be of some in- 
terest on very local scales). The entire process, in fact, is dominated 
by metal pollution and the strong yields of early, massive stars. 
We have to stress some caveats, though. Ejection of particles into 
the IGM is an unknown process. We have assumed winds origi- 
nated from stars (kinetic feedback), but different mechanisms (like 
gas stripping, shocks, thermal heating from stellar radiation, etc.) 
could play a role, as well, mostly at high redshifts, when objects 
are small and can easily loose part of their baryonic content. In 
addition, diffusion and conduction will probably alter the smooth- 
ness of metal and molecule distribution, but these phenomena have 



not been extensively studied, yet, and probably will depend on 
many parameters: e.g., the way metal or gas particles are ejected 
and mixed, how they are transferred away from the production 
sites, how strong is the efficiency of such processes, just to men- 
tion a few. Some attempts to address s uch issues have been done 
dSpitzerl|l962l: ICowie & McKeelll977l; lBrookshawl ll985l; ISarazin 
1988 1: iMonaghadl 19921; Iciearv & Monaghar 1999 ; lKlessen & Lir 
2003; Jubelgas et al. 2004; Monaghan et al. 2005; Wadsley et al 



2008; Greif et al, 2009), but much more realistic and detailed anal- 
yses are still strongly needed. 

To conclude, we have followed the structure formation process 
from very early times to first star formation and subsequent metal 
pollution. In the 1 Mpc side box simulations we have seen that, af- 
ter ~ 2 x 10 8 yr, molecular evolution leads the very first bursts 
of star formation (popIII), but metal enrichment is extremely fast 
(Figs.[T][2] andO in inducing the nowadays observed star formation 
mode (popII-I). In fact, we observe a steep increase of Z, with local 
values rapidly reaching and overtaking Z cr i t (Fig. [4}. Metal pollu- 
tion proceeds from the densest cores of star formation outwards, 
because of supernova ejections from high-density to lower-density 
environments (Figs. [2] and [3}. Rare, unpolluted regions can still 
survive, determining the simultaneous presence of the two star for- 
mation regimes, and Z cr u can affect the level of residual popIII star 
formation. As a result of this rapid pollution, we find that the aver- 
age contribution of the popIII component to the total star formation 
rate density is of a few times ~ 10~ 4 — 10~ 3 (with a maximum of 
~ 10 -2 — 10 _1 ) at 2 ~ 11 (Figs. [5|6| and|7l>. This general picture is 
preserved regardless of the precise value of the metallicity thresh- 
old Z cr it the slope of the popIII IMF, and their yields, but is quite 
sensitive to the popIII mass range. The change of the critical den- 
sity threshold for star formation influences little these conclusions, 
but can affect the onset and the decrement of the popIII contribu- 
tion to the total SFR (e.g. Fig. [§}. We have found similar results 
in larger-box simulations, either in the metal enrichment features 
(Figs.l9landll0t and in the star formation behaviour (Fig.ll It. 
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